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Abstract 



Background: Classical descriptions of enzyme kinetics ignore the physical nature of the in- 
tracellular environment. Main implicit assumptions behind such approaches are that reactions 
occur in compartment volumes which are large enough so that molecular discreteness can be ig- 
nored and that molecular transport occurs via diffusion. Though these conditions are frequently 
met in laboratory conditions, they are not characteristic of the intracellular environment, which 
is compartmentalized at the micron and submicron scales and in which active means of transport 
play a significant role. 

Results: Starting from a master equation description of enzyme reaction kinetics and assuming 
metabolic steady-state conditions, we derive novel mesoscopic rate equations which take into 
account (i) the intrinsic molecular noise due to the low copy number of molecules in intracel- 
lular compartments (ii) the physical nature of the substrate transport process, i.e. diffusion 
or vesicle-mediated transport. These equations replace the conventional macroscopic and de- 
terministic equations in the context of intracellular kinetics. The latter are recovered in the 
limit of infinite compartment volumes. We find that deviations from the predictions of classical 
kinetics are pronounced (hundreds of percent in the estimate for the reaction velocity) for en- 
zyme reactions occurring in compartments which are smaller than approximately 200nm, for the 
case of substrate transport to the compartment being mediated principally by vesicle or granule 
transport and in the presence of competitive enzyme inhibitors. 

Conclusions: The derived mesoscopic rate equations describe subcellular enzyme reaction ki- 
netics, taking into account, for the first time, the simultaneous influence of both intrinsic noise 
and the mode of transport. They clearly show the range of applicability of the conventional de- 
terministic equation models, namely intracellular conditions compatible with diffusive transport 
and simple enzyme mechanisms in several hundred nanometre-sized compartments. An active 
transport mechanism coupled with large intrinsic noise in enzyme concentrations is shown to lead 
to huge deviations from the predictions of deterministic models. This has implications for the 
common approach of modeling large intracellular reaction networks using ordinary differential 
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equations and also for the calculation of the effective dosage of competitive inhibitor drugs. 
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Background 

The inside of a cell is a highly complex environment. In the past two decades, detailed mea- 
surements of the chemical and biophysical properties of the cytoplasm have established that 
the conditions in which intracellular reactions occur are, by and large, very different than those 
typically maintained in laboratory conditions. One of the outstanding differences between in 
vivo and in vitro conditions, is that in the former, biochemical reactions typically occur in mi- 
nuscule reaction volumes [lj. For example, in eukaryotic cells, many biochemical pathways are 
sequestered within membrane-bound compartments, ranging from ~ 50nm diameter vesicles to 
the nucleus, which can be several microns in size [2]. It is also found that the total concentration 
of macromolecules inside both prokaryotic and eukaryotic cells is very large OS], of the order 
of 50 — 400 mg/ml which implies that between 5% and 40% of the total intracellular volume 
is physically occupied by these molecules [5] . The concentration of these crowding molecules is 
highly heterogeneous (see for example [6]), meaning that typically one will find small pockets 
of intracellular space, characterized by low macromolecular crowding, surrounded by a "sea" 
of high crowding; such pockets of space may serve as effective compartments where reactions 
may occur more easily than in the rest of the cytosol. Analysis of experimental data for the 
dependence of diffusion coefficients with molecular size suggests the length scale of such effective 
compartments is in the range 35-50nm [7], a size comparable to that of the smallest vesicles. The 
significant crowding also suggests that frequently an active means of transport such as vesicle- 
mediated transport, may be more desirable than simple diffusion as a means of intracellular 
transport. 

The volume of a spherical cavity of space of diameter 50nm is merely ~ 6.5 x 10 -20 liters, 
an extremely small number compared to the typical macroscopic reaction volumes of in vitro 
experiments (experimental attolitre biochemistry is still in its infancy - see for example [8]). 
These very small reaction volumes imply that at physiologically relevant concentrations (nano 
to millimolar), the copy number of a significant number of intracellular molecules is very small 
PQ and consequently that intrinsic noise cannot be ignored; for example 255fiM corresponds to 
an average of just 10 molecules in a 50nm vesicle and fluctuations about this mean of the order 
of 3 molecules [9J. 
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The traditional mathematical framework of physical chemistry ignores the basic physical 
properties of the intracellular environment. Kinetics are described by a set of coupled ordinary 
differential equations which implicitly assume (i) that the reaction compartment is so large 
that molecular discreteness can be ignored and that hence integer numbers of molecules per 
unit volume can be replaced by a continuous variable, the molar concentration. Since the 
number of molecules is assumed to be very large, stochastic fluctuations are deemed negligible 
and the equations are hence deterministic; (ii) the reaction compartment is well-stirred so that 
homogeneous conditions prevail throughout [9]. Both assumptions can be justified for reactions 
occurring in a constantly stirred reactor of macroscopic dimensions. However if diffusion is the 
dominant transport process inside the compartment then the homogeneity assumption holds 
only if the volume is small enough so that in the time between successive reactions, a molecule 
will diffuse a distance much larger than the size of the compartment. This comes at the expense 
of the first assumption. It hence appears natural that for intracellular applications, the first 
assumption, namely that of deterministic kinetics cannot be justified a priori. The second 
assumption can be justified if reactions are localized in sufficiently small parts of the cell and in 
particular for reaction-limited processes i.e. those for which the typical time for two molecules 
to meet each other via diffusion is much less than the typical time for them to react if they are 
in close proximity. For such conditions, a molecule will come within reaction range several times 
before participating in a successful reaction, in the process sampling the compartment many 
times which naturally leads to well-mixed conditions [TOl [TT] . 

In this article we seek to understand the magnitude of deviations from the classical ki- 
netic equations in small intracellular compartment volumes. We specifically focus on the case 
of reaction-limited enzyme reactions which allows us to relax the first assumption of physical 
chemistry while keeping the second one; this makes the mathematics tractable. We quantify 
deviations from classical kinetics in the context of the Michaelis-Menten (MM) equation; this is 
the cornerstone of present day enzyme kinetics and is a derivative of the traditional determinis- 
tic mathematical framework based on ordinary differential equations. In steady-state metabolic 
conditions, it is predicted to be exact. Thus this equation is ideal as a means to accurately test 
the effects of small-scale compartmentation on chemical kinetics. We consider three successive 
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biological models of intracellular enzyme kinetics, each one building on the biological detail and 
realism from the previous one (Figure 1). The models incorporate the intrinsic noisiness of 
kinetics in small compartments, the details of the substrate transport process to the compart- 
ment (diffusion or active transport) and the presence of intra-compartmental molecules other 
than substrate molecules which may modulate the enzyme-catalyzed process e.g. inhibitors. On 
the macroscopic level, i.e. for large volumes, the steady-state kinetics of all models conform 
to the MM equation. We test whether this equation holds on the on the length scale of small 
intracellular compartments by deriving the dependence of the ensemble averaged rate of product 
formation on the ensemble-averaged substrate concentration from the corresponding stochastic 
models in the steady-state. It is shown via both calculation and stochastic simulation that at 
these small length scales the MM equation breaks down, being replaced by a new more general 
equation. Practical consequences of this breakdown are illustrated in the context of the calcu- 
lation of the effective dosage of enzyme inhibitor drug needed to suppress intra-compartmental 
enzyme activity by a given amount. To make our approach accessible to readers not familiar 
with stochastic equations and their analysis, in the Results/Discussion sections we mainly focus 
on the biological/biophysical context and implications of the models together with the main 
mathematical results which are verified by simulation. Detailed mathematical derivations and 
the methods of simulation are relegated to the Methods section. 

Results 

Model I: Michaelis-Menten reaction occurring in a compartment volume of 
sub-micron dimensions. Substrate input into compartment occurs via a Pois- 
son process 

k 

This is the simplest, biologically-relevant case (Figure 1(A)). The reaction scheme is — ^ S + 

E v=^ C E + P. Substrate molecules (S) are continuously supplied inside the compartment 

fci 

at some rate k{ n , they reversibly bind to enzyme molecules (E) with rate constants ko (forward 
reaction) and k\ (backward reaction) to form transitory enzyme-substrate complex molecules 
(C) which then decay with rate ki into enzyme and product molecules (P). The substrate input 
is assumed to be governed by a Poisson process with mean h m \ this is consistent with substrate 
transport to the compartment being dominated by normal diffusion. The enzyme acts as a cat- 



6 



alyst, effectively speeding up the reaction by orders of magnitude. It is assumed that diffusion 
inside the compartment is normal and not rate-limiting on the catalytic process i.e. well-mixed 
conditions or rate- limited kinetics inside the compartment. Given these conditions we ask our- 
selves what is the relationship between the reaction velocity and the substrate concentration 
inside the compartment. The simplest approach consists of writing down the rate equations of 
traditional physical chemistry: 

[Br] = [E] + [C] = constant, (1) 
d[S]/dt = k in -k [E}[S]+k 1 [C], (2) 
d[C}/dt = k [E][S}-(ki+k 2 )[C], (3) 

^jf = UC]. (4) 

By imposing steady-state conditions we get the sought-after relationship which is simply the 
well-known MM equation: 

d[P] _ VmaxjS] , . 

dt K M + [sy {> 

where Km = (k\ + k2)/kn is the MM constant, v max = fef-Er] is the maximum reaction velocity 
and square brackets indicate the macroscopic concentrations. We note that steady-state condi- 
tions for substrate necessarily require that k{ n ^ Vmax otherwise the substrate will continuously 
accumulate with time. Though this approach is simple and straightforward, as mentioned in the 
introduction, the assumptions behind the formulation of the rate equations are not consistent 
with the known physical properties of the cytoplasm. In particular it is clear that if the volume 
of our compartment is very small (as is the case), the numbers of particles is also quite small, 
meaning that the concept of a continuous variable such as the average macroscopic concentration 
has little meaning. Rather we require a mathematical description in terms of discrete, integer 
numbers of particles and which is stochastic. The natural description of such cases is a master 
equation which is a differential equation in the joint probability function tt describing the system 

[si nana us]: 



^ = ^(©s 1 - 1)tt + ^(©^©c 1 " l)n S n E 7r 

+ hiQcQs 1 ~ l ) n CK + hiQcQp 1 - l)n C K, (6) 
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where ir = ir(ns, nc,np), ny is the integer number of molecules of type Y, £1 is the compartment 
volume, and O^ 1 are step operators defined in the Methods section. This equation cannot be 
solved exactly. However it can be solved perturbatively using the system-size expansion due to 
van Kampen |12j . This expansion is one in powers of the inverse square root of the compartment 
volume. In the Methods section, we calculate the first three terms of the expansion, namely 
those proportional to f2° and f2 -1 / 2 . The first term, being the dominant one for large 

volumes, gives back as expected, the rate equations Eqs. (l)-(4). The second term gives the 
magnitude of stochastic fluctuations about the macroscopic concentrations. Corrections to the 
rate equations and to the MM equation (due to small compartment volumes or equivalently 
due to intrinsic noise) are found by considering the third term. In the rest of the article, 
instead of using the reaction velocity v, we use the normalized reaction velocity, a, which is 
simply the velocity of the reaction, v, divided by the maximum reaction velocity, v max . Given 
some measured intracompartmental substrate concentration, [S*] = (n<j/f2) (angled brackets 
imply average), the relationship between the normalized reaction velocity predicted by the MM 
equation {olm = [&*]/ {Km + [S*])) and the actual normalized reaction velocity (a), as predicted 
by our theory, is given by: 



Hence the prediction of the MM equation is only correct, i.e. a = olm, in the limit of infinitely 
large compartment volumes, in which case the second term on the left hand side of Eq. (7) 
will become vanishingly small and can be neglected. For finite compartment volumes, the MM 
equation is not exact (except in the two limiting cases of olm ~ * and olm 1) but is at best 
an approximation, even though steady-state conditions are imposed; this is at odds with the 
prediction of the conventional deterministic theory. An inspection of Eqs. (7) and (8) shows that 
the magnitude of the deviations from the MM equation depends on the two non-dimensional 
quantities: (i) Km&, a measure of the rate at which enzyme-substrate combination events occur 
relative to the rate of decay of complex molecules; (ii) [-E^jfi, the total integer number of enzyme 
molecules in the compartment. 




(7) 



where 




(8) 
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As shown in the Methods section, the MM equation is found to implicity assume that the 
noise about the macroscopic substrate and enzyme concentrations is uncorrelated (this assump- 
tion has generally been found to be at the heart of many macroscopic models - for example 
see |16|): properly taking into account these non-zero correlations leads to the corrections en- 
capsulated by Eqs. (7) and (8). These correlations are expected to be small in two particular 
cases: (i) if Km is large; in this case when substrate molecules combine with an enzyme to 
form a complex, the latter dissociates very quickly back into free enzyme and thus successive 
enzyme-substrate events to the same enzyme molecule are bound to be almost independent of 
each other. The opposite situation of small Km would imply that the bottleneck in the catalytic 
process is the decay of complex rather than enzyme-substrate combination; if a successful com- 
bination occurs, the next substrate to arrive to the same enzyme molecule would have to wait 
until the complex decays, naturally leading to correlations between successive enzyme-substrate 
combination events, (ii) if the total number of enzyme molecules is large; in such a case, at any 
one time, the noise about the macroscopic concentrations will be the sum total from a large 
number of enzymes, each at a different stage in the catalytic process and each independent from 
all others, which naturally dilutes any temporal correlations. 

To estimate the magnitude of the deviations from the MM equation inside cells, we use 
the above two equations, Eqs. (7) and (8), to compute the absolute percentage error R p = 
100|1 — aM/ct\. These estimates are also computed by stochastic simulation of the master Eq. 
(6), using the exact stochastic simulation algorithm of Gillespie jlOj (see Methods for details 
regarding the method of simulation); this provides a direct test of the theory. Figure 2 shows 
the typical dependence of R p on c*m, as predicted by both theory (solid lines) and simulation 
(data points). Generally the agreement between the two is found to be very good; discrepancies 
increase as Km and compartment volume decrease but are small for parameter values realistic 
for intracellular conditions. The maxima of such plots gives the maximum absolute percentage 
error which is a measure of the maximum expected deviations from the MM equation. Table 1 
summarizes these estimates (theory and simulation) over wide ranges of the parameters typical 
of in vivo conditions: Km = lOfiM— 1000/iM [17], enzyme copy numbers of ten and one hundred 
per compartment which correspond to enzyme concentrations ranging from 4^M to 2.5mM and 
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compartment diameters ranging from 50nm to 200nm. Note that the maximum deviations from 
the MM equation are estimated to be less than approximately 20% and typically just a few 
percent over large ranges of parameter values - this robustness of the MM equation with respect 
to intrinsic molecular noise is indeed surprising, since strictly speaking it is only valid for infinite 
compartment volumes. 

The theory is always found to underestimate the actual deviations predicted by simulations; 
hence the theoretical expressions provide a quick, convenient way by which one can generally 
estimate a lower bound on the deviations to be expected from the MM equation without the 
need to perform extensive stochastic simulation. 

Model II: Michaelis-Menten reaction occurring in a compartment volume of 
sub-micron dimensions. Substrate is input into compartment in groups or 
bursts of M molecules at a time 

Model I captures the basics of a general enzyme-catalyzed process occurring in a small intracel- 
lular compartment. In this section we build upon this model to incorporate further biological 
realism. In particular, in the previous model we assumed that substrate input can be well de- 
scribed by a Poisson process, where one molecule at a time is fed into the compartment with 
some average rate ki n . This is the simplest possible assumption and approximates well the sit- 
uation in which molecules are brought to the compartment via normal diffusion. However there 
are many situations where this may not be the case; we now describe two such cases. 

The intracellular condition of macromolecular crowding limits the Brownian motion of molecules 
in the cytoplasm, this being reflected in the relatively small diffusion coefficients measured in 
vivo compared to those known in vitro for moderately to relatively large molecules. Experiments 
with inert tracer particles in the cytoplasm of Swiss 3T3 cells show that the in vivo diffusion 
coefficient is an order of magnitude less than that in vitro for molecules with hydrodynamic 
radius 14nm and diffusion becomes negligibly small for molecules larger than approximately 
25nm [7]; similar results have been obtained in Xenopus neurons |18j and skeletal muscle my- 
otubes [19J . If diffusion is considerably hindered, one expects active transport to become a more 
desirable mode of transport. Indeed there exists ample evidence for the active transport of 
macromolecules: they are typically packaged in a vesicle or a granule which is then transported 
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along microtubules or by some other means. It is also found that each vesicle or granule typ- 
ically contains several of these molecules (examples are: mRNA molecules - several estimated 
per granule [2Q|. I21j: cholesterol molecules which are transported in low-density lipoproteins [2] 
- approximately 1500 per lipoprotein). 

Generally an active means of transport is not exclusively linked with the transport of large 
substrate molecules. The cell being a highly compartmentalized and dynamic entity requires for 
its survival the precise transport of certain molecules from one compartment to another and a 
regulation of this transport depending on its current physiological state. Brownian motion leads 
to an isotropic movement of molecules down the concentration gradient and to a consequent 
damping of the substrate concentration with distance. In contrast active transport provides a 
directed (anisotropic) means of transport with little or no loss of substrate with distance, is 
independent of the concentration gradient and it is also easily amenable to modulation. 

Hence it follows that a more general process modeling molecular entry into an intracellular 
compartment is one in which M molecules are fed into the compartment at a rate k\ n \ the latter 
rate constant is the rate at which vesicles or granules arrive to the site of the compartment 
(Figure 1(B)). The total mean substrate input rate is then ki n = Mk® n . The special case 
M = 1 corresponds to Model I. We construct the relevant master equation and employ the 
system-size expansion as for the previous model (see Methods for details); it is found that the 
deterministic rate equations are exactly Eqs. (l)-(4) i.e. at the macroscopic level, given two 
reactions occurring in two different compartments, both with the same total mean substrate 
input rate k{ n but one occurring via diffusion (e.g. M = 1, k® n = 1) and the other via active 
transport (e.g. M = 10, fc^ = 0.1) , cannot be distinguished. However if the compartment 
volumes become small, then once again we find corrections to the MM equation and interestingly 
these corrections are sensitive to the mode of transport. The relationship between the normalized 
reaction velocity predicted by the MM equation {cum) and the actual normalized reaction velocity 
(a), as predicted by our theory, is given by Eq. (7) together with: 



This suggests that generally deviations from the predictions of the MM equation increase with 
the carrying capacity, M, of the vesicle or granule. To compare the effects of active transport 




(9) 
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and diffusion on the kinetics, we set M = 50 and adjusted k® n so that in all cases, the total mean 
substrate input rate for model II is equal to ki n , the input rate of Model I (i.e. the two models 
would be indistinguishable from a macroscopic point of view). Using the same procedure as for 
Model I, we computed the maximum percentage error using Eqs. (7) and (9) and also from 
simulations. The results are summarized in Table 2. Notice that now the deviations from the 
MM equation are much larger than before, running into hundreds of percent rather than the tens 
as for Model I. Because of the increase in substrate fluctuations, the quantitative accuracy of 
the theory is now less than before; it generally fares very well for compartments with diameters 
larger than ~ lOOnm and Km larger than ~ 100/iM. Nevertheless in all cases theory does 
correctly predict a large increase in discrepancy between the reaction velocities given by the 
deterministic MM equation and those from stochastic simulation compared to the case of Model 
I. The intuitive reason behind these increases in discrepancy is that substrate which is input in 
bursts enhances correlations between successive enzyme-substrate events. 

The explicit dependence of the reaction velocity on substrate concentration is complex and 
generally requires the solution of the cubic polynomial encapsulated by Eqs. (7) and (9). How- 
ever for small substrate concentrations, the equations simplify to a simple linear equation: 

a = [S*\ (k m 

Note that if the MM equation was correct, one would expect a = [S*]/Km- Indeed Eq. (10) 
reduces to the latter prediction in the limit of large volumes. Note also that this renormalization 
of the proportionality constant occurs only if the substrate input occurs in bursts, i.e. M > 1. 
These predictions of our theory are verified by simulations (Figure 3). 

Model III: Michaelis-Menten reaction with competitive inhibitor occurring 
in a compartment volume of sub-micron dimensions. Substrate input as in 
previous models 

In this last section, we further build on the previous two models by adding competitive inhibitors 
to the intracellular compartment in which enzymes are localized. A competitive inhibitor, /, 
is one which binds reversibly to the active site of the enzyme (forming a complex EI), thus 
preventing substrate molecules from binding to the enzyme and slowing down catalysis (Figure 
1(C)). In standard textbooks and in the literature, this is typically modeled by the set of reactions 
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M - 1 



2n(K M + [E T ])\ 



(10) 



(see for example [22]): S + E ^ C % E + P, E ^ EI, where h = h%[I\ and [I] is the 

fcl fc3 

inhibitor concentration. Note that it is implicitly assumed that inhibitor is in such abundance 
that the second-order bimolecular reaction between inhibitor and enzyme can be replaced by a 
pseudo first-order reaction with constant inhibitor concentration. We shall assume the same in 
our model. Substrate input into the compartment is considered to occur as in Model II since 
this encapsulates that of Model I as well. The deterministic model of this set of reactions leads 
to a MM equation of the form: 

d[ P ] _ VmaxjS] , , 

dt Km (i + p) + [sy 1 ] 

where (3 = [I]/Ki and Ki = k^/k^ is the dissociation constant of the inhibitor. The perturbative 
solution of the master equation describing the system is now significantly more involved than in 
previous models; the underlying reason for this is that the computation of the noise correlators 
to order Q° requires the inversion of a 6 x 6 matrix as opposed to a 3 x 3 one in previous models 
(see Methods for details). The analysis predicts corrections to the MM equation by postulating 
a new mesoscopic rate equation having the form of Eq. (7) together with: 

where q and di are coefficients with a complex dependence on the various enzyme parameters 
(these are given in full in the Methods Section). Table 3 shows the maximum percentage error 
computed using Eqs. (7) and (12) and also from simulations for the cases in which substrate 
input occur a molecule at a time and in bursts of 50 at a time. The parameter values chosen 
in the simulations and calculations (see caption of Table 3) are typical for many enzymatic 
processes: the bimolecular rate coefficients span the range 10 6 — 10 9 s _1 M _1 |22j . the backward 
decay processes are in the middle of the range 10 — 10 5 s _1 [22J, the inhibitor concentration is ten 
times larger than the total enzyme concentration (satisfying the implicit assumption that the 
inhibitor is in significantly larger concentration than free enzyme), and the intracompartmental 
enzyme concentrations are in the range 4 — 255/uM. The deviations from the MM equation in 
this case are more severe than in the previous two models, this being due to non-zero correlations 
between substrate and the complex EI in addition to the already present correlations between 
substrate and complex C. Note that the agreement between theory and simulations is overall 
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better than in previous models, even when the burst size is large, M = 50. As mentioned in the 
section for Model I, discrepancies between theory and simulation are generally found to decrease 
with increasing Km', for the case of competitive inhibition, the effective Km of the reaction 
is significantly larger than that of the enzyme (see Eq. (11)), which explains the increased 
agreement between theory and simulations for Model III compared to the previous two models. 

A significant number of drugs suppress a chain of biochemical reactions by reducing the 
activity of key enzymes in the pathway via competitive inhibition [T7j . The conventional method 
to estimate the required concentrations of these inhibitors involves plotting the variation of 
the enzyme activity with inhibitor concentration, [I], using the MM equation; the substrate 
concentration is kept fixed and is chosen so that at [I] = 0, the reaction velocity is close to 
the maximum, v max . Since there are significant corrections to the MM equation when reactions 
occur in intracellular compartments, it is not clear how accurate are estimates of [I] based 
upon it. Figure 4 compares the enzyme activity curve based on the MM equation with the 
theoretical predictions for the corrected enzyme activity curves based on the mesoscopic rate 
equation embodied by Eqs. (7) and (12), for compartments of diameter 50nm and lOOnm (inset) 
and for substrate input burst sizes of M = 20 and 50. The substrate concentration is chosen so 
that at [I] =0, v/vmax = 0.909 in all cases. We find that generally as the burst size increases, 
the actual inhibitor concentration needed to suppress enzyme activity by a given amount is 
larger than that estimated by the MM equation; this discrepancy decreases with increasing 
compartment volume. For the example in Figure 4, for the case in which substrate is input into 
the compartment in bursts of M = 50, the actual inhibitor concentration needed to decrease the 
enzyme activity from 0.909 to 0.1 is approximately 7 times larger than the MM estimate; if the 
compartment diameter is doubled (inset of Figure 4) , the actual inhibitor concentration needed 
is less than twice that of the MM estimate. Generally we find that for the typical parameter 
values of enzymatic reactions, the corrections to the enzyme-activity curves can be neglected for 
compartments larger than about 200nm in diameter. 
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Discussion and Conclusion 



In this last section we discuss some fine points regarding: (i)the assumptions behind the use of 
master equations which throws light on the range of use of the derived mesoscopic equations, 
(ii) the use of the system-size expansion to perturbatively solve the master equation and (hi) 
the assumption of steady-state metabolic conditions. We conclude by placing our work in the 
context of previous recent studies of stochastic enzyme kinetics and discuss possible experiments 
to verify some of the conclusions we have reached. 

We have implicitly assumed throughout the article that a single (global) master equation 
model suffices to capture the deviations from classical kinetics due to fluctuations in chemical 
concentrations inside a single subcellular compartment. As noted by Baras and Mansour [23], 
"the global master equation selects the very limited class of exceptionally large fluctuations 
that appear at the level of the entire system, disregarding important nonequilibrium features 
originated by local fluctuations." Hence the results presented here necessarily underestimate 
the possible deviations from classical kinetics, in particular the local fluctuations due to dif- 
fusion of molecules inside the compartment. These local fluctuations are typically small for 
reaction-limited processes (as in this article) but significant for diffusion-limited ones. To cap- 
ture them effectively, one would be required to spatially discretize the compartment into many 
small elements and describe the reaction-diffusion processes between these elements by means of 
a multivariate master equation [12 [ I23|. The latter is known as a reaction-diffusion master equa- 
tion; typically it does not allow detailed analytical investigation as for a global master equation 
and one is limited to stochastic simulation. Use of the global master equation is also restricted 
for compartments which are not too small: in particular the linear dimensions of the compart- 
ment should be larger than the average distance traveled by a molecule before undergoing a 
successful reaction with another molecule i.e. the length scale is much larger than that inherent 
in molecular dynamics simulation |23j . 

We have applied the systematic expansion due to van Kampen to perturbatively solve the 
master equation. It is sometimes a priori assumed that because this expansion is about the 
macroscopic concentrations, it cannot give information regarding the stochastic kinetics of few 
particle / small volume systems. This is true if one restricts oneself to the expansion to order 
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f2 i.e. the linear-noise approximation; this is commonly the case found in the literature since 
the algebra becomes tedious if one considers more terms. However we note that as argued and 
shown by van Kampen himself |12| . terms beyond the linear-noise approximation in the system- 
size expansion add terms to the fluctuations that are of order of a single particle relative to 
the macroscopic quantities and are essential to understanding how fluctuations are affected by 
the presence of non-linear terms in the macroscopic equation (substrate-enzyme binding in our 
case). In our theory we went beyond the linear-noise approximation. We find that the predicted 
theoretical results are in reasonable agreement, in many cases (comparison of bold and italic 
values in Tables 1, 2 and 3), with stochastic simulations of just a few tens of enzyme molecules 
in sub-micron compartments, which justifies our methodology. 

We have also imposed metabolic steady-state conditions inside the subcellular compartment. 
Technically this is convenient since in such a case one does not deal with complex transients. 
Also since under such conditions the MM equation is exact from a deterministic point of view, 
it provides a very useful reference point versus which to accurately compute deviations due to 
intrinsic noise. In reality one may not always have steady-state conditions inside cells, this 
depending strongly on the rate of substrate input relative to the maximum rate at which the 
enzyme can process substrate. Another possibility is that one is dealing with a batch reaction 
i.e. one in which a number of substrate molecules are transported at one go and just once to 
the subcellular compartment (e.g. via vesicle-mediated transport) and the reaction proceeds 
thereafter without any further substrate replenishment. This latter scenario is compatible with 
the presentation of the MM equation typical in standard physical chemistry textbooks. The 
MM equation is then an approximation (not exact as in steady-state case) to the deterministic 
kinetics, when substrate is present in much larger concentration than enzyme. This case is 
currently under investigation using the same perturbative framework used in this article. 

We note that this is not the first attempt to study stochastic enzyme kinetics. The bulk of 
recent studies \2A\ |25| 1261 [27] have focused on understanding the kinetics of a Michaelis-Menten 
type reaction catalyzed by a single enzyme molecule. Deviations from classical kinetics were 
found to be most pronounced when one takes into account substrate fluctuations [26] . These 
pioneering studies were restricted to a single-enzyme assisted reaction which reduces complexity 
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thereby making it ideal from a theoretical perspective; since the reaction is dependent on just 
a single enzyme molecule one also finds maximum deviations from deterministic kinetics. In 
reality, it is unlikely to find just one enzyme molecule inside a subcellular compartment - as 
mentioned in the introduction a physiological concentration of just a few hundred micromolar 
would correspond to few tens inside the typically smallest subcellular compartment. It is also the 
case that diffusion may not always be the main means of substrate transport to the compartment 
and that the reaction maybe more complex than the simple Michaelis-Menten type reaction of 
these previous studies. The present study fills in these gaps by using a systematic method 
to derive approximate and relatively simple analytic expressions for mesoscopic rate equations 
describing the kinetics of the general case of N enzyme molecules in a subcellular compartment 
with or without active transport of substrate and in the presence of enzyme inhibitors. Most 
importantly our approach shows the effects of intrinsic noise on the kinetics can be captured 
via effective ordinary differential equations. This enables quick estimation of the magnitude of 
stochastic effects on reaction kinetics and thus gives insight into whether a model or parts of 
a model should be designed to be stochastic or deterministic without the need for extensive 
stochastic simulation. In the present study, this approach enabled us to readily compute, for 
the first time, the deviations from deterministic kinetics for a broad range of realistic in vivo 
parameter constants (Tables 1, 2 and 3), a task which would be considerably lengthy if one had 
to rely solely on data obtained from ensemble-averaged stochastic simulations. 

We conclude by briefly discussing possible experiments which can verify the predictions 
made in this article. It is arguably not an easy task to perform the required experiments in real- 
time in a living cell. A viable alternative would consist of monitoring reaction kinetics inside 
single artificially-made vesicles. Pick et al [8] have shown that the addition of cytochalasin to 
mammalian cells induces them to extrude from their plasma membrane minuscule vesicles of 
attolitre volume with fully functional cell surface receptors and also retaining cytosolic proteins 
in their interior. The change in the intra- vesicular calcium ion concentration in response to 
surface ligand binding was measured using fluorescence confocal microscopy (FCM). Since the 
vesicle sizes are of typical small sub-cellular compartment dimensions (1 attolitre corresponds 
to a spherical vesicle of approximate diameter 120nm) and FCM allows the measurement of the 
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concentration of a fluorescent probe (via a calibration procedure), this experimental technique 
appears ideal to verify the predictions of Model I and of Model III for the case of diffusive 
substrate transport. Model II and Model III with vesicle-transport of substrate are probably 
much more challenging to verify since one then needs to construct the in vitro equivalent of 
microtubules. This is within the scope of synthetic biology and may be a possibility in the next 
few years. 

Methods 

We here provide full details of the calculations reported in the Results section. The system 
size-expansion which is at the heart of the analysis has to-date not been applied extensively to 
biological problems and thus we go into some detail in its elucidation in Subsection I, which 
is dedicated exclusively to Model I. For other recent applications of the general method in the 
context of reaction kinetics, see for example [28] and [29]. Subsections II and III (treating Model 
II and Model III, respectively) naturally build on the results of the first subsection and thus we 
only give the main steps of the calculations in these last two cases. Subsection IV has a brief 
discussion of the simulation methods used to verify the theoretical results. 

Model I: Michaelis-Menten reaction occurring in a compartment volume of 
sub-micron dimensions. Substrate input into compartment is modeled as a 
Poisson process 

The reaction scheme is S + E C E + P. The stochastic description of this system 

fci 

is encapsulated by the master equation which is a differential equation in the joint probability 
function tt describing the system: 

^ = h^O^ - 1)tt + ^(QsQc 1 ~ ±)nsn E ir (13) 
+ hiQcQ" 1 - l)n c vr + k 2 {Q c @ p l - l)n C TT, 

where tt = Tr(nc,np, ns), nx is the integer number of molecules of type X (where X = C, P, S), 
Q is the compartment volume, and B^ 1 are the step operators defined by their action on a 
general function g(nx) as: Q^ 1 g(nx) = g(nx ± 1). Note that the relevant variables are three, 
not four: the integer number of molecules of free enzyme (ng) is not an independent variable 
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due to the fact that the total amount of enzyme is conserved. The master equation cannot 
be solved exactly but it is possible to systematically approximate it by using an expansion in 
powers of the inverse square root of the volume of the compartments. This is generally called 
the system-size expansion |12j . 

The method proceeds as follows. The stochastic quantity, nx fluctuates about the macro- 
scopic concentrations [X]; these fluctuations are of the order of the square root of the number 
of particles: 

n x = n[X] + n 1/2 e x - (14) 

Note that since ng + nc = constant, it follows that nE = £l[E] — VL l l 2 ec- The joint distribution 
function and the operators can now be written as functions of the new variables, ex, giving: 
7T = U(ec,e P ,e s ,t) and G^ 1 = 1 ± Q-^d/dex + \®r l d 2 /de x + 0(fT 3 / 2 ); using these new 
variables the master equation Eq. (13) takes the form: 



dt 



1/2 fd[c] du , d[P] du , d[S] du\ _ nl/2 



+ 



+ 



dt dec dt dtp dt des 



- = nV2 ai n + o°a 2 n + n _1 / 2 a 3 n + o{^- v ) 



where 



(15) 



ai = -(k in + h[C] - k [E][S])— + ((h + k 2 )[C] - k [E][S})- k 2 [C] — , (16) 

oes oec oep 



id 2 .. 
a2= 2 km M + ~ A ^~ 



+ 



lf_d 

2\de s 
d 



_d_ 
de~ 



c 



(k [S][E] + kt[C}) + k 2 



d 



d 



d 



de s de c 



1 



[k (e s [E] - e c [S}) - hec] + -k 2 

2 \ dep dec 



dec dep 
d d 



e c 



[C\, 



(17) 



a3 



\(d_ d_ 

'2\d7^~de^ 
d 

k 



(koesiEj-koeciSj + hec) 



de s dec 



1/9 

esec + 2 k2 [-dT P 



d 



de 



c 



ec- 



(18) 



Note that in Eq. (18) terms which involve products of first and second-order derivatives, 
third-order derivatives or higher have been omitted - these do not affect the low-order moment 
equations which we will be calculating. 
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Analysis of fW 2 terms 

The terms of order fi 1 / 2 are the dominant ones in the limit of large volumes. By equating both 
terms of this order on the right and left hand sides of Eq. (15) and using Eq. (16), one gets the 
deterministic rate equations: 



d[S]/dt = k in -k Q [E][S]+k 1 [C], 
d[C}/dt = k [E}[S}-(k 1 + k 2 )[C], 
d[P]/dt = k 2 [C]. 



(19) 
(20) 
(21) 



These are exactly those which one would write down based on the classical approach whereby one 
ignores molecular discreteness and fluctuations. This is an important benchmark of the method 
since it shows that it gives the correct result in the limit of large volumes. On a more technical 
note, the cancelation of these two terms of order fi 1 / 2 makes Eq. (15) a proper expansion in 
powers of Q- 1 / 2 . By imposing steady-state conditions we have the Michaelis-Menten (MM) 
equation: 

d[P] _ v max [S] ^ 



dt K M + [S] ' 

where v max = fof-Er] is the maximum reaction velocity, [Er] = [E] + [C] is the total enzyme 
concentration which is a constant at all times and Km = (k\ + k2)/ko is the Michaelis-Menten 
constant. 

Analysis of 17° terms 

To this order, the master equation is a multivariate Fokker-Planck equation whose solution is 
Gaussian and thus fully determined by its first and second moments. The equations of motion 
for these moments can be straightforwardly obtained from the master equation to this order, 
leading to a set of coupled but solvable ordinary differential equations: 



(ec) 



-HE] 
k [E] 



ki + k [S] 
-k {K M + [S]) 



(es) 
(ec) 





r (4) i 




\ ^ 1 


dt 


(4) 


= A ■ 


(4) 




. (esec) _ 







(23) 



(24) 
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where, 

-2ko[E] 2(h + k [S] ) \ / k in + h [C] + k [S] [E] 

-2k (K M + [S}) 2k [E] \,B=[ k ([S][E] + K M [C}) | - (25) 

-2k 2 -2k 2 J V hn + k 2 [C] 

Note that the matrices and vectors in the above equations have been reduced to a simpler form 

by the application of a few row operations. Note also that these equations are independent 

of e v since the product-forming step is irreversible and hence the fluctuations in substrate and 

complex are necessarily decoupled from its fluctuations. At the steady-state, it is found that 

i e s,c) ~^ 0. From Eq. (14), it is clear that this implies that to this order the average number of 

substrate molecules per unit volume, (ris/Q), is simply equal to the macroscopic concentration, 

[S]. The same applies for complex molecules. Hence to this order in the system-size expansion 

there cannot be any corrections to the macroscopic equations or to the MM equation. By writing 

the macroscopic concentrations in Eqs. (24) and (25) in terms of ki n and solving, one obtains the 

variance and covariance of the fluctuations about the steady-state macroscopic concentrations. 

We here only give the result for the covariance since this will be central to our discussion later 

on: 

, N K M [E T }a 2 , v 

{eC€s) = K M + [E T ](l-ar (26) 

where a = ki n /v max is the normalized reaction velocity of the enzyme. 
Analysis of f2 -1 / 2 terms 

The system-size expansion is almost never carried out to this order because of the algebraic 
complexity typically involved, however it is crucial to find finite volume corrections to the de- 
terministic rate equations and in particular to the MM equation. Using the master equation to 
this order, the first moment of the complex concentration is governed by the equation of motion: 

d(e c )/dt = -k ([S] + K M ){e c ) + k [E](e s ) - k^ 1 ' 2 (e s e c ) . (27) 

Now the production of product P from complex occurs through a decay process which necessarily 
has to be described by a linear term of the form: ki n = ^(^c/^) (the steady-state condition). 
Since the steady-state macroscopic complex concentration is equal to [C] = ki n /k2, then it 
follows that to any order in the expansion we have (ec) = 0. This is always found to be the 
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case in simulations as well. Hence it immediately follows from Eq. (27) that the average of 
fluctuations about the macroscopic substrate concentration are non-zero and given by: 

{€s) = aUWY ( } 

From a physical point of view, this indicates that the steady-state concentration of substrate in 
the compartment is not equal to the value predicted by the MM equation (i.e. [S]) and hence 
the non-zero value of the average of the fluctuations about [S] . The real substrate concentration 
inside the compartment is obtained by substituting Eqs. (28) and (26) in Equation (14), leading 
to: 

/I^\ = rQi + KmGi2 (OQ\ 

\ n I L J ^ (1 - a) [K M + [E T }(1 - a) 2 ]n { ' 

An alternative mesoscopic rate equation replacing the MM equation 

The renormalization of the steady-state substrate concentration indicates the breakdown of the 
MM equation; this phenomenon occurs because of non-zero correlations between noise in the 
substrate and enzyme concentrations, (e^ec), which the MM equation implicity neglects. To 
obtain the alternative to the latter, one needs to obtain a relationship between the normalized 
reaction velocity, a and the real substrate concentration {ns/ty; writing [S] in terms of a and 
substituting in Eq. (29), one obtains this new relation: 

a+(l- {Hs/n) ) f(a)n-i - {ns/Q) (30) 
+ V K M + (n s /n)) na) ' 1 -K M + (n s /ny m 

a 2 

= K M + [E T]{ l-a? (31) 

Note that in the limit of large volumes, the second term on the left hand side of Eq. (30) becomes 
vanishingly small and we are left with the MM equation. In the results section the quantity on 
the right hand side of Eq. (30) is referred to as au since this is the normalized reaction velocity 
which would be predicted by the MM equation given the measured substrate concentration 
{ns/ty inside the compartment. A quick estimate of the magnitude of the error that one 
is bound to incur by using the conventional MM equation can be obtained by substituting 
a = 1/2 (i.e. enzyme is half saturated with substrate) in Eqs. (30) and (31), solving for oim and 
then using this value to compute the fractional error e = 1 — au/oi. This leads to the simple 
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expression: 

e = [l + n([E T ]+4K M )}- 1 (32) 

We finish this section by noting that Eq. (30) will be found to be valid generally and not only for 
the simple Michaelis-Menten scheme treated in this section; the details of the reaction network 
come in through the form of Eq. (31) which is reaction-specific. 

Model II: Michaelis-Menten reaction occurring in a compartment volume of 
sub-micron dimensions. Substrate is input into compartment in groups or 
bursts of M molecules at a time. 

A natural generalization of Model I which has direct biological application is when substrate 
molecules are fed into the compartment M at a time with mean rate k® n . The total mean 
substrate input rate is then equal to ki n = Mkf n . The master equation for this process is Eq. 
(13) with the first term on the right hand side replaced by Q(Og M — l)k® n . This leads to the 
following change in the expression for 02 (Eq. 17): 

1 d 2 1 d 2 

TT^-m 77~9 > ^ki n M 9 . (33) 

2 deg 2 de z s 

Note that since the expression for a\ (Eq. 16) is unchanged, the deterministic equations are 
precisely the same as those of Model I. However now the fluctuations about the macroscopic 
substrate concentration are enhanced by a factor M; consequently the entries in the vector B 
in Eq. (25) need the change ki n — > ki n M. The analysis proceeds in the same manner as before. 
The mesoscopic rate equation replacing the MM equation is now given by Eq. (30) together 
with: 

a[a + ±(M-l)] 
K M + [-EtJ(1 - a) 



The fractional error rate evaluated at a = 1/2 gives: 

M 



(35) 



M + n([E T ] + AK M ) 

This clearly shows that generally larger deviations from the predictions of the MM equation are 
expected in this case compared to those computed for Model I. 
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Model III: Michaelis-Menten reaction with competitive inhibitor occurring in 
a compartment volume of sub-micron dimensions. Substrate input as in two 
previous models. 

Competitive inhibition is modeled by the set of reactions: — ^ S+E C — ► E+P, E EI, 

ki fc3 

where k^ = k®[T\ and [I] is the inhibitor concentration (similar models have been studied by 
Roussel and collaborators [30|, [3Tj in the context of biochemical oscillators though these assume 
M = 1). In the rest of this section, we change the notation of enzyme-inhibitor complex from EI 
to V, just to make the math notation easier to read. The substrate input into the compartment 
is considered to occur as in Model II since this encapsulates that of Model I as well. The master 
equation for this system is: 

^ = kln(Q s M - 1)tt + ^(©se^ 1 - l)n s n E n + ^(©c©^ ~ l ) n c* 

+ k 2 {Q c Q p l - l)n c ir + k 3 (O v - l)n v + k^Oy 1 - l)n E . (36) 

The change of variables from nx to ex is done as before, however note that now the conservation 
law for enzyme is different than in the two previous models. The total enzyme concentration 
is now equal to [E T ] = [E] + [C] + [V] from which it follows that n E = Q,[E]- f2 1/2 (e c + e v ). 
The description is chosen to be in terms of numbers of molecules of types C, S and V and thus 
E being a dependent variable does not show up explicitly in the step operators of the master 
equation above. 

Due to the significant number of changes in the terms of the expansion from those of previous 
models, we will show the equivalent of Eqs. (15)-(18) in full. The master equation in the new 
variables ex is given by: 



m nl/2 /d\c] &u_ d\p] &n_ d\s] an d[v) an 



dt " \ dt dec dt dep dt des dt dey, 

o^Vn + n°a 2 u + n~ 1 / 2 a 3 n + oin- 1 ) (37) 



where 



01 = - (ki n + h[C] - k [E)[S])— + ((h + k 2 )[C] - k [E)[S])— (38) 

oes dec 

+ {k,[V}-UE])4--k 2 [C)- d 



dey dep'' 
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(39) 



If d d 
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d d 



de s de c 



i, / s ^ V i, d 2 i, a 2 , 

+ 2 /C2 fe"^J ec+ 2" 3 94^-2" 4 94 (ec + ev ' ) - (40) 
Analysis of f] 1 / 2 terms 

As for previous models, these terms give the macroscopic equations. Equating both terms of 
this order on the right and left hand sides of Eq. (37) and using Eq. (38), one obtains: 

d[S] fdt = k in - k [E] [S] + h [C] , (41) 

d[C]/dt = k [E][S] - (h + k 2 )[C], (42) 

d[P]/dt = k 2 [C], (43) 

d[V]/dt = k 4 [E] - k 3 [V]. (44) 

In the steady-state we have the Michaelis-Menten (MM) equation: 

d[P] = VmaAS] , . 

dt K M {l + l3) + [Sy 1 ' 

where (3 = [I\/Ki and Ki = k^/k® is the dissociation constant of the inhibitor. 

Analysis of 17° and SI -1 / 2 terms 

The equations for the first moments are easily obtained and we shall not reproduce them here; 
suffice to say that at steady-state, it is found that (es,c,v) ~^ which implies that to this order 
in the system-size expansion there cannot be any corrections to the macroscopic equations or 
to the MM equation. The addition of a new species, V, does however substantially increase the 
algebraic complexity in the equations of motion for the second moments computed using terms 
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up to order 0°. In particular the matrix A is now a 6 x 6 matrix, rather than the 3x3 matrix 
of the previous two models. 



dt 



r <4> i 
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(47) 



and 



5 = 



(48) 



/ kinM + hiCl + koiSHE] \ 
k ([S][E] + K M [C}) 
h[E] + k 3 [V] 



±(fc in M + fc 2 [C]) 

In the above equations we have defined k! = k^ + k^. Note also that the system of equations has 
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been simplified through the application of a few row operations. 

Now to next order, i.e. O -1 / 2 , the first moments of the concentrations of molecules of type 
C and V are governed by the equation of motions: 

d(e c )/dt = -k ([S] + K M )(e c ) - k [S](e v ) + k Q [E](e s ) - k Q- l / 2 {{e s e c ) + (e s e v )) (49) 
d(e v )/dt = -k 3 (e v ) - k 4 ((e c ) + (e v )) (50) 



As in previous models, since the production of product P from complex occurs through a decay 
process, it follows that at steady-state, (ec) = which also implies (ey) = from Eq. (50). 
Hence it follows from Eq. (49) that (es) = [{es^c) + {esev)]/^ l ^[E\- The two cross correlators 
can be estimated to order Q° by solving Eqs. (46)-(48). The non-zero value of (es) implies a 
renormalization of the substrate concentration inside the compartment and hence to a new rate 
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equation replacing the MM equation. This is obtained exactly in the same manner as previously 
shown for Model I. The mesoscopic rate equation is found to be given by Eq. (30) together with: 

where the numerator coefficients are given by: 

c o = + k 3 (0 + l) 3 K 2 M k o [E T ], (52) 

ci = + K M ((3 + l) 2 [(/3 + l)[E T ]kt - (3/3 + 2)[E T ]k K M k 3 + k o 0v max K M ], (53) 

c 2 =- K M ((3 + 1) [2(2/3 + 2 + l)[E T ]k 2 3 - ((30 2 + 4(3+ l)[E T ]k K M + (54) 

- 0(0 + l)v max - k [E T ] 2 )k 3 + /3(1 + 2(3)k v max K M -(J3+ l)[E T }k v max }, 

c 3 = + [(1 + 3/3 + /3 3 + 2>p 2 )[E T )K M ki - (/3(/3 + l) 2 [E T ]k K 2 M + (55) 
(-0(13 + lfv max - 2(1 + (3)[E T } 2 k )K M + /3(/3 + l)[E T ]v max )k 3 + 
/3 2 (1 + (3)k v max K 2 M - (2 + 3/3 + 2p 2 )[E T ]k v max K M }, 

c 4 = - [(-(0 + l)[E T ] 2 k K M + f3((3 + l)[E T ]v max )k 3 + (56) 

- [E T }((3 + (3 2 + l)k v max K M }, 

and the denominator coefficients by: 

d = + K 2 M k k 3 (l + f3)\ (57) 

d 1 = + K M k 3 (f3 + lf[0(k 3 - k K M ) + k 3 ], (58) 

d 2 = + K M k o (0 + l) 2 [k 3 [E T }(0 + 2)+ v max ], (59) 

d 3 = + (/3 + l)[k 2 2 [E T ] - kop 2 k 3 K M [E T } + 2k 2 f3[E T ] (60) 

- k o k 3 0K M [E T ] - k (3v max K M + k 2 [E T ]}, (61) 
d 4 = + [E T ]k [k 3 [E T ] + k 3 (3[E T ] + 

Vmax 

]. (62) 

Note that Ylt=o c « = such that at a = 0, there is no correction to the MM equation i.e. 
ajw = also. The case = reduces to Model II, i.e. /(a) is given by Eq. (34). 

Stochastic simulation 

In this section we briefly describe the simulation methods used to verify the theoretical results 
which are described in detail in the Results section. All simulations were carried out using 
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Gillespie's exact stochastic simulation algorithm, conveniently implemented in the standard 
simulation platform, Dizzy |32j . 

The data points in Figure 2 were generated by iterating the following four-step procedure: 
(i) pick a value for a between and 1. This gives the substrate input rate ki n = av max ; (ii) run 
the simulation and measure the ensemble-averaged substrate concentration, (n s /0) = [S*] at 
long times; (hi) compute «m using the MM equation, oim = [S*]/{[S*] + Km)', (iv) compute the 
absolute percentage error R p = 100|(1 — aj\f/a)\. The solid curves in Figure 2 were obtained by 
numerically solving the cubic polynomial in a given by Eqs. (7) and (8) in the Results section 
for given values of an and then using the above expression for R p . Figure 3 is generated in the 
same manner as Figure 2, except that: in step (i) we fix M and pick a value for a between 
and 1. Since k% n = Mk® n , the required simulation parameter is k® n = av max /M; step (iv) is not 
needed. The solid curves were obtained by numerically solving the cubic polynomial in a given 
by Eqs. (7) and (9) in the Results section for given values of [S*]. The y-axis for this figure 
is v/v max = om for the MM equation and v/v max = a for the stochastic model. Figure 4 is 
obtained by numerically solving the quintic polynomial in a given by Eqs. (7) and (12) in the 
Results section together with the coefficients given by Eqs. (52)-(62) in the present section; the 
inhibitor concentration, [I], is varied while the substrate concentration, [S*], is kept fixed. The 
substrate concentration is chosen so that at [I] = 0, v/v max = 0.909 in all cases. Note that for 
models I and II, a M = [S*]/([S*] + K M ) while for Model III, a M = [S*}/([S*] + (1 + P)K M ). 
Note that the error bars are very small on the scale of the figures and thus are not shown. 
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Table 1: Maximum Percentage error in reaction velocity from prediction of the MM equation for 
Model I. The copy number indicates the total number of enzyme molecules per compartment. 
Values in bold and in square brackets are those estimated by simulation; the italic values are 
obtained from the derived theoretical expressions, Eqs. (7) and (8). 



D/nm 


K M = 10/xM 


100/xM 


lOOO^M 


Copy No. 


50 


11.83 [17.00] 


4.09 [4.33] 


0.59 


10 


100 


4.74 [5.00] 


0.73 [0.74] 


0.08 


10 


200 


0.90 


0.10 


0.01 


10 


50 


3.98 [5.33] 


1.88 [2.02] 


0.43 


100 


100 


2.10 [2.23] 


0.52 [0.52] 


0.07 


100 


200 


0.61 


0.09 


0.01 


100 



Table 2: Maximum Percentage error in reaction velocity from prediction of the MM equation for 
Model II. The copy number indicates the total number of enzyme molecules per compartment. 
Values in bold and in square brackets are those estimated by simulation; the italic values are 
obtained from the theoretical expressions, Eqs. (7) and (9). 



D/nm 


K M = 10/xM 


100//M 


lOOO/iM 


Copy No. 


50 


225.40 


152.83 [291.56] 


45.43 


10 


100 


161.59 [331.66] 


52.74 [58.39] 


6.82 [6.99] 


10 


200 


65.09 


8.45 [8.50] 


0.88 


10 


50 


32.97 


30.17 [61.66] 


18. 14 


100 


100 


30.78 [66.03] 


19.76 [24.52] 


5.57 [6.06] 


100 


200 


21.27 


6.61 [6.91] 


0.85 


100 
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Figure 1: Schematic illustrating the three models considered in this article. (A) Model 
I: Michaelis-Menten reaction occurring in a compartment volume of sub-micron dimensions 
(shown by dashed rectangle). Substrate input into compartment occurs via a Poisson process 
i.e. diffusion-mediated substrate transport. (B) Model II: As for Model I but now substrate 
is input into compartment in groups or bursts of M molecules at a time i.e. vesicle-mediated 
substrate transport along microtubules (MT). (C) Model III: Michaelis-Menten reaction with 
competitive inhibitor (7) occurring in a small subcellular compartment. Substrate transport as 
in previous two models. 
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Figure 2: Deviations from the predictions of the MM equation for diffusion-mediated 
substrate transport. (Model I) Plot of the Percentage Error in reaction velocity, R p = 
100 1 1 — oim/oi\, versus the normalized reaction velocity of the MM equation, olm for 10 enzymes 
(green) and 100 enzymes (red) with Km = 10/iM in compartments with diameter lOOnm (A) 
and 50nm (B) . The solid lines show the theoretical predictions, as encapsulated by Eqs. (7) 
and (8); the data points are obtained by stochastic simulation (see Methods for details). 
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Figure 3: Deviations from the predictions of the MM equation for vesicle-mediated 
substrate transport. (Model II) Testing the validity of the MM relationship at small substrate 
concentrations for the case in which substrate input into compartments occurs in bursts. The 
data is for 10 enzymes with Km = 100/xM in compartments of diameter (A) 200nm (circles), 
(B) lOOnm (diamonds) and (C) 50nm (crosses); substrate is input M = 50 molecules at a time. 
The deterministic prediction for all three cases is the same MM equation shown by the green 
curve. In contrast, the stochastic models, [Eqs. (7) and(9)], predict different rate equations for 
each case (red solid lines). Data points are obtained by stochastic simulation (see Methods for 
details). Note that v/v max = oim and a for solid green and red lines respectively. 
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Figure 4: Effects of intrinsic noise on the inhibition of enzyme activity in small 
compartments. (Model III) Plots of normalized enzyme activity versus normalized inhibitor 
concentration (measured in units of the total enzyme concentration [Er]) for 10 enzymes with 
Km = 100/iM in compartments of diameter 50nm and lOOnm (inset). The colors correspond 
to: (red) MM equation; (green) stochastic model, M = 20; (blue) stochastic model, M = 50. 
The latter two curves are those predicted by theory [Eqs. (7) and(12)]. Parameters same as 
mentioned in caption of Table 3 (except for [I] , which is a variable in the present case) . Substrate 
concentrations chosen so that at [I] = 0, v/v max = 0.909 in all cases. Black dashed lines contrast 
the inhibitor concentration required to decrease enzyme activity from 0.909 to 0.1 as predicted 
by the MM equation and the stochastic models. Note that v/v max = au and a for solid red 
and blue/green lines respectively. 
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Table 3: Maximum Percentage error in reaction velocity from prediction of the MM equation for 
Model III. The total number of enzyme molecules per compartment is ten in all cases. Values 
in bold and in square brackets are those estimated by simulation; the italic values are obtained 
from the theoretical expressions, Eqs. (7) and (12). The parameters are: ko = 10 9 s _1 M _1 , k\ = 
k 3 = 1000s- 1 , kl = 10 7 s- 1 M- 1 , and [I] = 10[E T ]. 



D/nm 


K M = 10/xM 


WOfiM 


lOOO/iM 


M (burst size) 


50 


67.8 [76.8] 


67.8 [76.5] 


67.8 


1 


100 


20.8 [26.4] 


20.6 [26.1] 


20.6 


1 


200 


2.8 


2.7 


2.7 


1 


50 


1001.8 


234.9 [169.4] 


86 


50 


100 


343.7 [345.5] 


73.4 [75.2] 


26.2 [31.5] 


50 


200 


71.4 


11.3 [11.5] 


3.6 


50 
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